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Abstract. The problem of Laplacian growth is considered within the Loewner-equation 
framework. A new method of deriving the Loewner equation for a large class of growth problems 
in the half-plane is presented. The method is based on the Schwarz-Christoffel transformation 
between the so-called 'mathematical planes' at two infinitesimally separated times. Our method 
not only reproduces the correct Loewner evolution for the case of slit-like fingers but also can 
be extended to treat more general growth problems. In particular, the Loewner equation for 
^ ' the case of a bubble growing into the half-plane is presented. 
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1. Introduction 



The Loewner equation pQ is an important result in the theory of univalent functions [2] that 
has found important applications in nonlinear dynamics, statistical physics, and conformal field 
theory [31 0]. In its most basic formulation, the Loewner equation is a first-order differential 
equation for the conformal mapping gt(z) from a given 'physical domain,' consisting of a 
complex region P minus a curve I\ emanating from its boundary, onto a 'mathematical domain' 
represented by P itself. Usually, P is either the upper half-plane or the exterior of the unit 
circle, but recently the Loewner equation for the channel geometry was also considered [5]. The 
Loewner equation depends on a driving function, here called a(t), that is the image of the growing 
tip under the mapping gt{z). An important development on the theory of the Loewner equation 
was the discovery by Schramm [6] that when the driving function a{t) is a Brownian motion the 
resulting Loewner evolution describes the scaling limit of certain statistical mechanics models. 
This result spurred great interest in the so-called stochastic Loewner equation [3]. 

Recently, the deterministic Loewner equation was also used to study the problem of Laplacian 
fingered growth in both the half-plane and radial geometries [8] as well as in the channel 
geometry jo] . In this case, the driving function a(t) has to follow a specific time evolution in 
order to ensure that the finger tip grows along the gradient lines of the corresponding Laplacian 
field. The idea of using iterated conformal maps to generate aggregates was first deployed by 
Hastings and Levitov [9j in the context of stochastic growth models, such as diffusion limited 
aggregation; see, e.g., [10] for further developments along those lines. A deterministic version 
of the Hastings-Levitov model that is closely related to the Loewner-equation approach — albeit 
not using explicitly such a formalism — was studied in [TT] . 

In this paper we consider the problem of Laplacian growth within the context of the Loewner 
evolution, and present a new method of deriving the corresponding Loewner equation for a broad 



class of growth models in the half-plane. Our method is based on the Schwarz-Christoffel (SC) 
transformation between the mathematical planes gt and gt+ T , where r is an infinitesimal time 
interval. More specifically, the method consists of expanding the integrand of the SC formula in 
powers of the appropriate infinitesimal quantity (related to r) and then performing the integrals 
up to the leading-order term. Our method correctly yields the Loewner evolution for the case 
of slit-like fingers studied before [5j [8] . More importantly, the method is able to handle more 
general growth problems, so long as the growth rule can be specified (in the mathematical plane) 
in terms of a polygonal curve, in which case the Schwarz-Christoffel transformation can be used. 
An example is given for the case of a bubble growing from the real axis into the upper half-plane. 

2. Loewner Equation for Slit-like Fingers 

2. 1 . The Case of a Single Finger 

In order to set the stage for the remainder of the paper, we wish to begin our discussion by 
considering the simplest Loewner evolution, namely, that in which a curve starts from the real 
axis at t = and then grows into the upper half-z-plane H, where 

H = {z = x + iy € C : y > 0}. 

The curve at time t is denoted by T t and its growing tip is labeled by j(t). Now let gt(z) be 
the conformal mapping that maps the 'physical domain,' corresponding to the upper half-z- 
plane minus the curve Tt, onto the upper half-plane of an auxiliary complex w-plane, called the 
'mathematical plane,' i.e., we have w = gt(z), where 

g t : m\T t -> H, (1) 

with the curve tip j(t) being mapped to a point a(t) on the real axis in the w-plane; see Fig. [TJ 
Furthermore, we consider the growth process to be such that the accrued portion of the curve 
from t to t + r, where r is an infinitesimal time interval, is mapped under gt{z) to a vertical 
slit in the mathematical w-plane; see Fig. [TJ The mapping function gt{z) must also satisfy the 
initial condition 

9o{z) = z, (2) 

since we start with an empty upper half-plane. We also impose the so-called hydrodynamic 
normalization condition at infinity: 

g t (z) =z + 0(- n ), z^oo. (3) 
\z\ 

These conditions specify uniquely the mapping function gt{z). 

From a more physical viewpoint, the problem formulated above belong to the class of 
Laplacian growth models where an interface evolves between two phases driven by a scalar 
field (ft(x,y;t), representing, for example, temperature, pressure, or concentration, depending on 
the physical problem at hand. In one phase, initially occupying the entire upper half-plane, the 
scalar field (ft satisfies the Laplace equation 

V 2 = 0, (4) 

whereas in the other phase one considers 0=const., say (ft = 0, with the curve Tt representing 
a finger- like advancing interface between the two phases. (Here the finger is assumed to 
be infinitesimally thin.) The complex potential for the problem can then be defined as 
w(z,t) = ift(x,y;t) + i(ft(x,y;t), where tft is the function harmonically conjugated to (ft. On 
the boundary of the physical domain, consisting here of the real axis together with the curve Tt, 
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Figure 1. The physical z-plane and the mathematical w- and ("-planes at times t and t + r, 
respectively, for a single finger in the upper half-plane. The mapping gt maps the curve V (at 
time t) onto a segment of the real axis on the it;-plane, whereas the accrued portion of the 
curve (during time interval r) is mapped to a vertical slit. The mapping gt is obtained as the 
composition of gt+ T and the slit mapping F; see text. 



we impose the condition = 0, whereas at infinity we assume a uniform gradient field, V0 ~ y, 
or alternatively, 

w(z, t) pa z, z — > oo. (5) 

From this point of view, the mapping function gt{z) introduced above corresponds precisely to 
the complex potential w(z, t) of the problem. In particular, the fact that in the w-plane the curve 
grows along a vertical line implies that the finger tip grows along gradient lines in the z-plane. 
To specify completely a given physical model, one has also to prescribe the interface velocity, 
which is usually taken to be proportional to some power 77 of the gradient field: v ~ | \70 1 77 . For 
most of the problems considered here the specific velocity model is not relevant, in the sense 
that the finger shapes will be independent of the exponent r/, which only affects the time scale of 
the problem. (However, there are situations, such as the case of competing asymmetrical fingers 
[5], where different 77's may yield different patterns.) 

For convenience of notation, we shall represent the mathematical plane at time t + r as the 
complex (-plane and so we write ( = gt+ T (z). Now consider the mapping w = F((), from the 
upper half-(-plane onto the mathematical domain in the w-plane; see Fig. [H The mapping 
function gt+ T (z) can then be given in terms of gt(z) as 

gt+T = F~ 1 og t , (6) 

where F -1 is the inverse of F(Q. The above relation governs the time evolution of the function 
gt(z) and naturally leads to the Loewner equation. A standard way of showing this is to construct 
the slit mapping F(Q explicitly, substitute its inverse in (|6|), and then take the limit r — > 0. 
One then finds the Loewner equation 

9t(z) = . f l) m , (7) 
g t {z) - a(t) 

where d(t) is the so-called growth factor which is related to the tip velocity. One can show [5] 
that d(t) = |//'(a(t))| _r '/ 2 ~ 1 , where ft{z) is the inverse of gt(z). Here, however, the specific form 
of d(t) is not relevant since we can always rescale the time coordinate in (|7|) so as to set d(t) = 1. 
From symmetry, one also gets 

d(t) = 0, (8) 



so that a(t) = ao = const., which implies that the tip j(t) simply traces out a vertical line in 
the z-plane. 

In more general situations, such as the case of multiple fingers, the function -F(C) can no 
longer be obtained in closed form and so one has to resort to alternative approaches to derive 
the Loewner evolution. Previous methods [U E] consider a series of compositions of the basic 
one-slit mapping. Here, however, we will apply a more direct method based on the Schwarz- 
Christoffel transformation to obtain the Loewner equation for slit-like fingers in the half-plane. 
In the next section, our method will be extended to treat the case of a growing bubble in the 
upper half-plane. To illustrate how the method works, let us first use it to derive the Loewner 
equation ([7]) for a single finger. 

We begin by inverting ([6]) so as to write 

g t = Fog t+T . (9) 

The mapping -F(C) from the upper half-^-plane onto the upper half-w-plane with a vertical slit 
(see Fig. [T]) is easily found by a direct application of the Schwarz-Christoffel formula [12]. One 
then finds 

f 9t+T C — a(t + T~\ 
g t = F(g t+T ) = / y=k 1 + =dC + a(t), (10) 

where Co = a(t + t) + v. Note from Fig. [I] that the parameter v is related to the (infinitesimal) 
height of the slit in the w-plane. The above integral can be performed exactly, as already 
mentioned, but here we take an alternative approach, namely, we first expand the integrand in 
(fTUj) in powers of v and then compute the relevant integrals afterwards. To do that we first 
rewrite f|10[) in the form 

{■9t+ T Ar 
9t= u 577 T. + «(*)■ ("■) 



, Co vi-^ 2 [c-^ + t)]- 2 

After expanding the integrand in powers of v 2 and performing the relevant integrals, one obtains 
that, up to order v 2 , equation (flT]) becomes 

v 2 l2 

9t+r-9t = '. , v +a(t + r)-a(t). (12) 

g t+T - a{t + r) 

Now expanding this equation up to the first order in r, dividing by r, and then taking r — > 0, 
one gets 

d(t) , . 

9t = , + ' m + a, 13 
g t (z) - a(t) 

where 

d(t) = \\mj-. (14) 

t^O IT 

Using the boundary condition lim 9t _j. 00 g t = 0, which follows from ([3]), yields precisely the 
Loewner equation (J7|) together with the condition (JS]). 



2.2. The Multifinger Case 

Here we consider the case of multiples fingers Vi, i = 1, n, growing from the real axis into the 
upper half-z-plane, as shown in Fig. [2l As before, the map gt(z) maps the physical domain in 
the z-plane onto the upper half-ui-plane and the tips are required to grow along gradient lines, 
so that the accrued portions of the curves Tj during an infinitesimal time r are mapped under 
gt(z) to vertical slits emanating from the real axis; see Fig[2l 
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Figure 2. The physical and mathematical planes for the case of multiple fingers in the upper 
half-plane. 

The mapping w = F((), from the upper half-(-plane to the upper half-w-plane with n vertical 
slits, can again be easily obtained from the Schwarz-Christoffel transformation |12j : 



9t = F(g t+T ) = [J C ai(t + r) = ^ + ^ ^ 

where Co = a j{t + r) + Uj for a given j. [We remark parenthetically that in writing (|15p we have 
assumed, for simplicity, that the slits in the tu-plane are mapped under F~ 1 (w) onto symmetrical 
segments on the real (-axis; see Fig [21 Rigorously speaking, this is valid only in the limit that 
the slit heights become vanishingly small, that is, when r — > 0, which is the relevant limit for us 
here.] 

The integral in (|15p cannot be performed exactly for arbitrary n, hence in order to obtain 
the Loewner equation for this case we first need to expand the integrand in powers of the 
infinitesimal parameters Vi and then proceed with the integration. Note, however, that each 
term in (|15p is of the same form as that appearing in (|10p for the case of a single finger. We 
can thus build upon our experience with that case to treat the present situation. In particular, 
we notice that the mixed terms involving different fj's in the expansion of the integrand in (|15p 
are all of orders higher than v 2 and hence need not be considered, for they do not contribute to 
the final result in the limit r — > 0. Thus, to the extent that the mixed terms can be neglected, 
we can rewrite (1151) as 



n rgt+T 



i=l^° Jl-V?[C-Oi(t + T)]* 



Now repeating the exactly same procedure used for the single finger case, see Eqs. (Ilip - (ll3p . 
one readily obtains 

« di( t) _ - d (t) 

t^9t-ai(t) a 3 {t) - ai {t) 3Kh K J 

where 



v} 



di{t) = f T (18) 



After using the condition limg t _s. 0O 

gt = in ([T7]), we get the Loewner equation for multiple 

curves 

9t = jZ^%-y ( 19 ) 
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Figure 3. The physical and mathematical planes for the case of a growing bubble in the upper 
half-plane. Here the new interface at time t + r is mapped under gt to a tent-like shape; see text 
for details. 



with the time evolution of the points ai{t) being given by the following system of ordinary 
differential equations 

ai{t) - aj(t) 

If the growth factors di(t) are all the same, we can again rescale the time variable so as to 
set di = 1. In particular, in the case of two symmetrical fingers (i.e., d\ = c?2 = 1), equation 
(]19p can be integrated exactly to yield the mapping function gt(z), from which the finger shapes 
can be computed analytically [5j. A related exact solution for two fingers was obtained in |13j . 
An alternative derivation of (I19p was given elsewhere [8] using a composition of n single-slit 
mappings. Our method, however, is somewhat more direct in the sense that it considers a single 
mapping with n slits as shown in Fig[2l In the next section we will extend our method to include 
the case of a growing bubble in the half-plane. 



3. Loewner Equation for a Bubble in the Half-Plane 

Here we consider the problem of an interface starting initially from a segment, say, [—1, 1], along 
the real axis and then growing into the upper half-z-plane, as indicated in Fig. [3j As before, we 
denote the interface at time t by Tj and consider the mapping gt(z) from the physical domain 
in the z-plane to the mathematical w-plane: 

g t : H\r t -> H. 

The mapping function gt(z) must also satisfy the hydrodynamic normalization condition ([3]). 
The growth dynamics is specified by requiring that each point on the interface grow along 
gradient lines, but with a growth rate that is maximum at the center (i.e., at the bubble tip) 
and zero at both 'contact points,' z = ±1, in such a way that the interface at time t+T is mapped 
under gt[z) to a tent-like shape, as shown in Fig. [3l The height of this tent, or alternatively, 
the angle a formed with the real axis, is linearly related to r in the limit r — > 0, as we will see 
below. Here we consider only the case of symmetric bubbles, where the contact points, z = ±1, 
are mapped by gt(z) to the points w = +a(t), respectively, and the tip (point B in Fig. [3]) is 
mapped to the origin. 

Since the domain in the tf-plane has a polygonal shape, the mapping w = -F(C) can once 
again be obtained from the Schwarz-Christoffel transformation. In this case we have 

rgt+r /-2a /-k 

* = F ^' ) = U } K>-k + W- « + *>- (21) 
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Figure 4. Loewner evolution for a bubble growing in the upper half-plane for d(t) = 1. The 
bubble starts at t = from the interval [—1, 1] on the real axis. The different curves represent 
the interface at time intervals of At = 0.3. 

which up to the first order in the parameter a becomes 



Ja<t+r) ( V J 

After performing the integral above and taking the limit r — > 0, one readily obtains the following 
generalized Loewner equation: 



d(t) = 1 by rescaling the time coordinate, in which case we have a(t) = oo4*. Recall also that 
the Loewner equation (|23p must be supplemented with the initial condition go(z) = z. 

In Fig. 2] we show numerical solutions of f|23|) with ao = 1 for various final times i, starting 
from t = 0.3 up i = 2.1 with a time separation of At = 0.3 between successive curves. To 
generate the curves shown in Fig. H] we used the numerical scheme described in |13] . More 
specifically, we start with a 'terminal condition' gt = w, for w £ [— a(t), a(t)], and integrate 
the Loewner equation ([23]) backwards in time, using a Runge-Kutta method of second order, to 
get the corresponding point go on the interface. (See also |14J for a recent review on numerical 
integration of the Loewner equation.) From Fig. [J] one sees that the bubble initially grows 
somewhat slowly and then rapidly expands and tends to occupy the whole plane for large times. 
In fact, one can show that (for d = 1) the tip velocity grows exponentially with time as t — > oo. 
It is possible to modify the growth factor d(t) so to have the tip velocity related to the gradient 
of the field 4>(x,y), as discussed in Sec. but this does not change the interface shapes and only 
alters the time scale of the bubble evolution. 




(22) 



g t {z) = d(t) {[g t + o(t)] ln[g t + a(t)} + [g t - a{t)} ln[g t - a(t)} - 2g t lng t } 
together with the governing equation for the logarithmic singularity a(t) 

a(t) = (In 4)d(t)a(t), 
where the growth factor d(t) is now given by 



(23) 



(24) 




4. Conclusions 

We have presented a novel method to derive the Loewner equation for Laplacian growth 
problems. The method is based on the Schwarz-Christoffel (SC) transformation and consists 
of expanding the integrand of the SC formula in the appropriate infinitesimal parameter, 
performing the relevant integrals, and then taking the limit in which the infinitesimal parameter 
goes to zero. Our method is able to reproduce the Loewner evolution for the problem of slit-like 
fingers in both the half-plane (Sec. [2]) and the channel geometry (not shown here). Furthermore, 
the method can be extended to treat more complicated growth problems, so long as the growth 
dynamics in the complex-potential plane can be specified in terms of a polygonal boundary, 
in which case the Schwarz-Christoffel transformation can be used. We note however that the 
requirement that the growth rule be formulated in terms of a polygonal curve is not as restrictive 
as it seems, for any simple curve can in principle be approximated by a piecewise linear function. 
[Such more general growth models are currently under investigation.] In particular, we obtained 
the Loewner equation for a novel situation in which a bubble grows from a segment of the 
real line into the upper half-plane. Although in this case we refer to the evolving interface as 
a growing 'bubble,' in contrast to the slit-like 'fingers' of Sec. [2j this terminology should not 
be taken literally. Depending on the physical problem at hand, such growing interface may 
represent, say, an expanding front in combustion experiments or in electrochemical deposition 
[5]. Of course, further work is necessary to relate more directly the growth models discussed 
here to experiments. 
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